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Abstract 

Using the Resonant Valence Bond (RVB) concept and the Schwinger-boson mean field approx- 
imation, we investigate a two dimensional anisotropic Heisenberg antiferromagnet. We find that 
increasing the coupling ratio (a = J y /J x ), a disordered-ordered phase transition appears beyond 
a critical value (a c ). When a > a c there exists a finite XJv beyond which an ordered-disordered 
phase transition occurs, while for a < a c there always exists a gap. Also we find that a decou- 
pling temperature Tp exists for ID and 2D isotropic case, when T approaches Tjj the RVB order 
parameter decreases rapidly to zero. 
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I. INTRODUCTION 



Low dimensional Heisenberg antiferromagnet (HAF) system has always been an interest- 
ing problem to theoretical physicists. Many new concepts and techniques are developed in 
this field. In the 1950's P.W.Anderson and Kubo developed the spin-wave theory (SWT) jlj|. 
Spin-wave theory predicted for the existence of an ordered ground state and gaplessed spin- 
wave excitations above it for higher dimensional magnets. In 1973 2], P.W.Anderson pro- 
posed his well known Resonant Valence Bond (RVB) model to describe some spin-1/2 sys- 
tems. Then in 1983 [jj, Haldane proposed his conjecture, it pronounced that the excitation 
spectrum of a linear-chain Heisenberg antiferromagnet (LCHA) with integer spin has a finite 
energy gap Eh above its singlet ground state, while a LCHA with half-odd-integer spin has 
a gapless spectrum. Later Affleck, Kennedy, Lieb and Tasaki (AKLT) 4| generalized the 
RVB model to spin-1 LCHA by proposing the Valence-Bond-Solid (VBS) concept to treat 
the conjecture. They said that these valence bonds are formed by two 1/2 spins as a sin- 
glet tl — IT on the nearest neighbors, while the two 1/2 spins on the same site should be 
symmetrized to form a triplet state S — 1. Since then many theoretical and experimental 
works are made and now the existence of the gap in integer spin systems is generally 

accepted. 

Most recently the two dimensional (2D) anisotropic Heisenberg antiferromagnet prob- 
lem has stired many people's interests. Many theoretical technics (spin-wave theory g 
Schwinger-boson mean field theory |9j and nonlinear a model |lfj ]) are used to study this 
problem. It seems that there exists complex phase transition in these anisotropic systems due 
to the coupling ratio (a- = J y /J x ) and temperatue (/3 = ksT). In this article, we will apply 
the Schwinger boson mean field theory (SBMFT) to treat the problem. And a compariation 
to other methods will be made. In Section II, we present the self consistent equations and 
make a brief discussion. In Section III, we calculate the self consistent equations at ground 
state numerically and find a disordered-ordered phase transition in the spin-anisotropy phase 
diagram. In Section IV, we discuss the Neel transition (AF ordered-disordered phase) due to 
the temperature, and a comparison to exponents fl Q on CMC, wil, be made. Also, 
we find that a decoupling temperature To exists for ID and 2D isotropic case, when T > Td 
the valence bond will be decoupled, and the self consistent equations have no solution in 
this region. 
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II. DERIVATION OF SELF-CONSISTENT EQUATIONS 



A. SBMFT Derivation 

The 2D square lattice Heisenberg antiferromagnet can be given by the following Hamil- 
tonian j^, 

H = J x Si ■ Sj + J y Si ■ S m , (1) 

<i,j> <l,m> 

where the sums < i,j > and < l,m > are defined for the nearest neighbors along the x 
and the y directions. J x > is the exchange coupling along the x direction and J y > along 
the y direction. Without losing the generality, we suppose the coupling ratio a = J y / J x is 
not larger than 1. If we take the limit a <C 1, the problem degenerates into the ID case. 
When a increases, we can expect a ID to 2D crossover. 

We employ the Schwinger-boson mean field theory (SBMFT) developed by Arovas and 



Auerbach to treat the Hamiltonian(l). Unlike the Holstein-Primakoff transformation in 
the spin-wave theory, we use the Schwinger-bosons to represent the spin operator. 
With the Schwingger-boson operators: s-f = a, s± = b, the spin S can be written as: 

S+ = a + b, S~ = b + a (2a) 
S z = \ (a + a - b + b) (2b) 
S = \ (a + a + b + b) . (2c) 

It's easy to establish that these operators fulfill the restriction: [S z , S + ] = S + , 

[s z ,s-] = -s-. 

Since the relation(2), we can rewrite the Hamiltonian(l) as: 

H = J X J2 Si-Sj + Jy S l -S m + i^J2{ a t a i + K b i- 2S ) ( 3 ) 

<i,j> <l,m> q 

Where the index q runs over the whole lattice, a^a q and b+b q are the number operators 
of spin-up and spin-down bosons. And /i is the lagrange multiplier to conform that there are 
25* Schwinger-bosons per site. Since the spin values are all the same, the lagrange multipliers 
for different sites are equal. The Hamiltonian(3) has the translational symmetry and we can 
apply the Fourier transformation (FT) to it. 



To simplify the discussion , let's study the x direction only, the H x term can be written 

as: 



H X = J X Y^ (Si ■ Sj - \n inj ) + J, \ n * n i + 2 ( a t a i + K h i - 2S ) ( 4 ) 
<*>i> <i,j> q I — I 

Now let's introduce the RVB order parameter A (or the valence bond for spin-1 case) |l2j |: 
A z = \ (bia i+z — CLib i+z ). A z is the order parameter along z(x or y)-direction, and z is the 
unit vector along z(x or y)-direction. Then the Hamiltonian has the following form: 



H x = C - 2 (f + J X S) NS - J x { «xbf - bt +x at) A x + h.c.} + J x ^ (f + J X S) n t , 

<i,j> i 

(5) 

where C = 2J X N \A X \ 2 + J X S 2 N, N: number of lattice sites. To simplify the derivation 
we ignore the C term in Hamiltonian(5), for it does not involve \x. It only shifts the whole 
system's energy with a certain value, and does not influence the excitation's spectrum. Thus, 
the Hamiltonian(5) can be rewritten in momentum space as: 

H x = ^2{\ x (a+a fe + bt k b-k) + Ikx^btk + l* kx a k b-k} - 2NS\ X (6) 

k 

Hence the Hamiltonian(3) can be expressed as: 

H = H X + H y = ^{\ (a+a fe + b± h b- k ) + 7*4^ + - 2NSX, (7) 

k 

with: X x = w+Y* 8 ) ; \ y = (/J+2 2 JyS) , A = X x + X y ; ^ kx = 2iJ x A x sin k x , >y ky = 
2iJ y A y sin k y , 7 fc = ^ kx + ^ ky . 

Using the following Bogliubov transformation [Tj 



h = -(^t) V2 e^at + ( *£>) V2 e^(3 k , (8b) 



the Hamiltonian can be diagonalized as: 

Hmf = J2 i Ek ( a t a * + \)+E k (PtP k + l)}-(25+ 1) iVA (9) 
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with the energy spectrum: E k = ^X 2 — 7^7^. 

Using the diagonalized Hamiltonian(9), we obtain the free energy: 



F = -{3- 1 In {Tr exp {-(3H)} = | ^ In {2 sinh } - (2S + 1) NX. (10) 

k 

Thus we get the self consistent equations for A x and A y : 

A * = Y1 ( coth ^r) {Wf) ( J - A - sin + sin *») ( lla ) 

k 

A y = J2 (coth^ffc) (f^f) (J^sin^ + J^sinAg. (lib) 



From the condition g = §f = (n| or = 25 iV |l2j, we obtain the third self- 

i 

consistent equation: 2£±1 = ^ (E^N)' 1 (coth ^pO- And these self-consistent equations can 
be written as integrations: 



A x = y (f) 2 ( c °th^) (^) (J x A x smk x + JyAySinky) (12a) 

A y = / (f) 2 ( coth/ ir) {^) (J x A x smk x + JyAySinky) (12b) 

^ = / (f) 2 (^)"Mh^). (12c) 

where the self consistent equations integrate from — ir to n, and Ek = 
\J X 2 — 4 ( J X A X sin k x + J y A y sin fc^) 2 . 



B. Spin-spin correlation 

Since Si • Sj can be expanded as: 

Si ■ Sj = -\ (aj~bf - bjaf) (bidj - dibj) + \ (afaj - bfbj) (a^a { - b+bi) , (13) 

where Rj = Ri + R, R = R x + R y . Using expression (13), the spin-spin correlation can be 
expressed as jll| : 

(Si ■ Sj) = (S ■ S R ) = \f (R)\ 2 - \g (R)\ 2 , (14) 
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TABLE I: renormalized factor Fr for differrnt spin values: the SBMFT solution of large spin 
system (such as spin-1) is more believable than the small spin system (spin-1). 



spin 


1/2 


1 


S 


CO 


Fr 


3/4 


8/9 


AS(S + 1)7(25 + l) 2 


1 



with: 



f(R) 



(9 ) e -k- coth ^-y^ 



^(^) = l/(f) 2 ^tcoth^. 



(15a) 
(15b) 



Consider R = 0, (Si ■ Sj) = (Sf) = S (S + 1) (h — 1), we can renormalize the spin-spin 
correlation function as: (So • Sr) = Fr (\f (R)\ 2 — \g (-R)| 2 ), Fr is the renormalized factor. 

From equation(14) it's easy to conclude that if (i,j) belongs to the same sublattice g(R) 
vanishes, and if (i,j) belongs to the different sub attice f(R) vanishes. That means there 
exists antiferromagnet (AF) order in this system ^, 
both f(R) and q (R) vanish. This indicates the AF order is short ranged and dereases 
expotentialy 0, Q . 
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III. GROUND STATE PROPERTIES 



A. ID solution 



First, we shall study the a = ID limit and a = 1 2D isotropic limit. In the ID case 
the self consistent eqs. can be written as: 



. dk coth ^ J X A X sm 2 k x 
A x = / 2 * 16a 

2S+1 fdk coth Sf. 

— = ./ (16b) 



2 

x) • 



with: Ek = y A 2 — 4 ( J x A x sin k 
Consider the ground state T — > 0, i.e., coth — > 1. We can calculate A x and A 
numerically for different spin values. 
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TABLE II: ID SBMFT solution for different spin values at ground state: Spin-0 case only calculated 
as a reference, and itself has no specific physical meaning. Both spin-1/2 and spin-1 showed an 
energy gap (2S factor in Eh formula indicates that there are 2S bosons excited for a physical 
excitations). 



spin 


A x . 


A 


Eh = 2Suj m in 








0.5J X 





1/2 


0.6792 


1.38 J x 


0.2A25J X 


1 


1.1816 


2.3647J X 


0.1701.4 


oo 


oo 


2 J 3u 






Experiments showed that only integer spin case has a finite energy gap (the Haldane 
gap). This implies that the SBMFT does not fit for the ID half-odd-integer spin problems. 
Theoretically this is because the topological term in the long-wavelength effective action. 
A formalism based on fermionic representation for the spin operators is relevant for 
the half integer case. Numerical result for spin-1 case is agree with other authors' result. 
Experimently the gap (0.41J.,;) [7J is larger than the SBMFT result, this is because our 
model is too simple. The model for a real spin-1 system contains a single- ion anisotropy 

i 

B. 2D isotropic solution 

Similarly we can study the self-consistent eqs. of 2D isotropic case (l^ . We assume 
Jx = Jy = J and obtain the 2D isotropic self consistent eqs.: 



, f dk\ 2 coth ^f^JA sin fc x (sink x + sinfc„) 
A. = A„ = A= / l^-j > (17a) 

2S+1 { (dk\ 2 coth^Si 



with: Ek = \J X 2 — 4 J 2 A 2 (sin k x + sin k y ) 2 . For the ground state T — >• and coth — > 
1, there exists an finite maximum value of 25" + 1 for the 2D integration(17). So there exists 
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an critical spin value S c : 



{2S + 1) 



max 



2S C + 1 



/ (i) 2 (v^) 



(18) 



If S < S c , there exists a gap, and if S > S c gapless, a Bose-Einstein condensation occurs at 
the zero energy point Q,^^. Numerical calculation shows that S c = 0.1956. This means 
that 2D isotropic Heisenberg antiferromagnets on a square lattice with nearest coupling 
always have ordered ground states, and gapless. D. Yoshioka has studied the problem 
for a spin-1/2 isotropic square lattice using the slave-fermion mean field theory. Likewise, 
we can generate to an anisotropic case. Since the spectrum is gapless, the Bogoliubov 
transformation(8) lost its meaning at certain ki in which E(ki) = 0. D. Yoshioka introduces 
a special transformation to diagonalize the mean field Hamiltonian(7) at Ki. 




(19a) 




(19b) 



and Cki, Cfci fulfill the following commutate relations: 




[U.£]=°- 



(20b) 



(20a) 



Hence we obtain: 




(21a) 




(21b) 



[( ki ,H}= [(+,H} =0. 



(21c) 
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Now the self consistent eqs. should be rewritten as: 



N ^ 



a coth 



A, 



N ^ 



sin k x (J X A X sin k x + J y A y sin k y 
A 



E k 



A coth £§t 



'1 - 4,/tJ + no5jfc,fej 



1 ^-^ sin fc^ ( J X A X sin fc x . + Jj,A y sin k y) 



N 



A 



Acoth^fk. . 

^ (1 - kM ) + n d kikt 



and the spin-spin correlation for the gaplessed spectrum: 



(22a) 
(22b) 
(22c) 



(S -S R ) = F R (\f(R)\ 2 -\g(R)\ 2 ) 



a coth 

- (1 - Ok,ki) + n 6 k ,k x 



yK ' IN ^ A 



Eh 



A coth £ft 



1 - 4,fcJ + n 8 kiki 



(23a) 
(23b) 

(23c) 



Compared eqs. (23) with eqs. (15), we can conclude that if there is no gap then there is 
Neel long range order, and if there is a gap there only exists short range order and long 
range disordered. 



C. 2D anisotropic solution 

As we have discussed in section III. A there exist significant difference between ID and 
2D isotropic Heisenberg antiferromagnet systems. There should exist a critical coupling a c , 
when a < a c , the system behaves like a ID chain (quasi ID spin system) and when a > a c 
it has a Neel order and_gaplessed spectrum (2D spin system). This prediction has been 
made by T.Sakai et.al. |8( using the spin wave theory and M.Azzour |9| using the SBMFT. 



Recently A.Parola et.al. [13] declare that for spin-1/2 system there exists a disordered 
transition induced by anisotropy at about a < 0.1, they find that the disordered phase is 
gapless and its wavelength can be interpreted in terms of a decoupled ID chains. 

In M.Azzour et.al.'s article 9], they provide a solution of SBMFT for large spin approx- 
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TABLE III: critical coupling a c for different spin values at ground state T = 0. 





This article 


T.Sakai [8] 


M.Azzour [9J 


D.Senechal [1£] 


Method 


SBMFT 


SWT 


SBMFT 


NLcj model 


spin- 1/2 


0.13 


0.03367 


0.02 


0.55 


spin-1 


0.009 


0.0013 


0.0019 


0.03 



imation: 



And D.Senechal et.al. 



r« -2ttS 

a c = be 



work out a nonlinear a model solution: 



(24a) 
(24b) 



a,. 



4exp 



87T 2 

3(at) 2 



(25) 



where they choose at = 2.1. 

Employing the condition Eh = and T — > 0. We discuss the ordered-disordered phase 
transition of ground state by calculating the self consistent eqs.(12). Compared with other 
authors' results, our numerical results are listed in table III. The ordered-disordered phase 
diagram for different spin value and anisotropy is given in Fig. I. 



IV. EXCITED STATE PROPERTIES 
A. The Neel phase transition 

Experimentally ID spin systems are realized in quasi-lD compounds in which the coupling 
J x is much higher than the coupling J y . It is the coupling ratio a = J y /J x determines the 
degree of "quasi-one-dimensionality" of the materials. For instance, a is estimated to be 
0.02 in CsNiCk, and NENP [Ni(C 2 H 8 N2)2N0 2 ( y ClO^} 0.0006. And experiments indicate 
that there exists notable difference between them. In CsNiCU ID behavior, the existence of 
gap is 0b se„ed above a c riti ca, Mature (the Nee, te m pe ratUr e T N ) of aW 5K QQ. 
While in NENP ID behavior, the existence of gap is observed at temperature as low as can 
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be reached. This indicates that for spin-1 system a c is not larger than 0.02. The result 
from nonlinear a model in table III seems too large [lOj, while SBMFT solution and SWT 
solution are permitted. 

Calculating the self consistent eqs.(12) numerically, we can find the T/v for CsNiCl%. 
It approximates IK, much lower than the experimental result. Perhaps it is due to the 
compound C sNid^s structure is ABX3 type, not cubic jf L And for a real system, the 
Hamiltonian should include the anisotropic term D (S?) 2 a, Q D is produced by the 
coupling of a spin to the anisotropic orbital motion. 

r~l 

R.Botet et.al. [14] found that a gap exists for 1.6 > y > —0.5, and the gap de- 
creased rapidly for negative D. Experiments indicate that the anisotropy constant is negative 
j- = —0.038 5]. So the calculated T/v (D = 0) should be smaller than the T/v for a real 
system (D = —0.038 J). The ordered-disordered phase diagram for different spin values and 
temperature with a = 0.02 is given in Fig. 2. 



B. The temperature-anisotropy phase diagram of spin-1 system 

Numerical solution of eqs.(12) indicate that when (3 —>■ 0.46 the RVB order parameter A 
decreases rapidly to zero for both ID and 2D isotropic spin-1 systems. Suppose A x = A y = 0, 
we can calculate the critical temperature Tp, denoting that when T > Tp the boson pairs 
on the nearest neighbors will be decoupled. 

coth ^- = 2S + 1 (26) 

Solving the equation (26), we find that To = 2.16 for spin-1 case of a = 1 and a = 0. 
But, for the anisotropic case < a < 1, the self consistent eqs.(12) has no solution under 
the condition = A y = or A x 7^ 0, A y = 0. This indicates that when increasing the 
temperature, the Valence-Bond-Solid concept lost its validity before the boson pairs on the 
nearest sites decoupled. The phase diagram of spin-1 Heisenberg antiferromagnet can be 
divided into three regions: the ordered phase, the disordered phase, and the decoupled phase. 
Valence-Bond-Solid concept lost its validity in the decoupled phase, we must consider other 
order parameters. The phase diagram for spin-1 2D anisotropic Heisenberg antiferromagnet 
is given in Fig. 3. 
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For spin- 1/2 systems such as KCuF 3 and Cu-0 layers etc., T D is calculated out to be 
0.910. 

V. SUMMARY 

Using the the Schwinger-boson mean field theory, we have investigated the phase transi- 
tion in anisotropic 2D Heisenberg antiferromagnet for various spin values. We find that the 
anisotropic HAF system can be characterized as two types according to different coupling 
ratio a = J y /J x . When a < a c , such as NENP for spin-1 system, it belongs to the quasi 
ID spin systems. For the quasi ID spin systems, the disordered phase always exists. When 
a > a c , such as CsNiCl^ for spin-1 system, it belongs to the usual 2D spin system. For 
the usual 2D spin system, the disordered phase only exists when T > T/y. So we can use 
CsNiCl 3 and other ABX 3 type compounds (such as RbNiCl 3 , CsNiF 3 et.al.) for testing 
anisotropic 2D HAF problem and use NENP for testing Haldane conjecture and other ID 
spin-chain HAF problems. Also we find that for the ID and 2D isotropic HAF system there 
exists another critical temperature T D , when T > T D , A x — A y — 0, the boson pairs will 
be decoupled. But for the anisotropic case (0 < a < 1), there has not such a decoupled 
solution for A x = A y = or A x ^ 0, A y = 0. This implies order parameters other than 
RVB should be considered in the high temperature region. 
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FIG. 1: Spin value - anisotropy a phase diagram at ground state T = 0. 
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FIG. 2: Spin value - temperature Jx/ksT phase diagram with a = 0.02. 
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FIG. 3: Temperature Jx/ksT - anisotropy a phase diagram for spin-1 2D anisotropic Heisenberg 
antifer r omagnet . 
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